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Definitions of fractional derivatives as fractional powers of derivative operators are 
! suggested. The Taylor series and Fourier series are used to define fractional power of 

\Q ■ self-adjoint derivative operator. The Fourier integrals and Weyl quantization proce- 

dure are applied to derive the definition of fractional derivative operator. Fractional 
generalization of concept of stability is considered. 
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1 Introduction 

7— I ■ 
>. 

The theory of integrals and derivatives of non-integer order goes back to Leibniz, Liouville, 
Riemann, Grunwald, and Letnikov [HE]- Fractional analysis has found many applications 
in recent studies in mechanics and physics. The interest in fractional equations has been 
growing continually during the last few years because of numerous applications. In a short 
period of time the list of applications becomes long. For example, it includes the chaotic 
£3 dynamics OH], material sciences [5J, [61 [TJ, [8] , mechanics of fractal and complex media [9[ [TUl 

[TT] , quantum mechanics [121 [32] , physical kinetics [3J EH E51 EES] , plasma physics [TO [20] , 
electromagnetic theory [TS], [TS], [20], astrophysics [21], long-range dissipation [221 [20], non- 
Hamiltonian mechanics [231 [2H 12H] , long-range interaction [28l [29l [30] . 

It is known that we can define a fractional power of operator [331 EH ESI EH [321 EH EZ] . 
The integer power of operator can be easily realized. Therefore, we can realize the fractional 
power as the integer power series. In this paper, we use equations that represent the fractional 
power as a series of integer powers series. This representation allows us to define the fractional 
power of operator as a series of integer powers of operator. As a result, we obtain the 
definition of fractional derivatives as a fractional power of derivative operator. 

Note that the well-known Riemann- Liouville fractional derivative can be represented as 
a power series of derivatives of integer order pQ: 

°° d n 
^a + = $^n0r,a,Q0 — , 

n=0 
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where 

{-l^aTin - a) {x - a) n ' a 
n(X ' a ' a) ~ r(l-a)r(n + l) T(n + l-a) 

for the functions that are analytical in the interval (a, b). 

In Sec. 2, we point out some well-known definitions of functions of bounded and un- 
bounded operators. In Sec. 3, the fractional derivatives are defined as fractional powers of 
coordinates that considered as Taylor series. In Sec. 4, the fractional derivatives are con- 
sidered as fractional powers of coordinates that described as Fourier series for the interval. 
In Sec. 5, the fractional derivatives are defined as fractional powers of coordinates by using 
the Fourier integrals. In Sec. 6, the fractional derivatives are defined as fractional powers 
of coordinates by using the Weyl quantization. In Sec. 7, using fractional derivatives, we 
define the stability with respect to fractional variations. 

2 Function of bounded and unbounded operators 

Let us point out some well-known definitions of functions of bounded and unbounded oper- 
ators [SDESlEalElEH]. 

2.1 Power series 

Let us consider a bounded linear operator A that is defined on the linear space E, and 
A G L(E,E), where L(E,E) is a space of linear maps of E. Suppose the function f(x) is 
an analytical function of the variable x such that it can be represented as a power series 

oo 

f(x) = J2fnX n . 
n=0 

Then, we can define 

oo 

f(A) = J2fnA n . 
n=0 

The operator f(A) is a linear bounded operator A on space E 
function of operator is defined by 

oo 1 

e A = Y -A n . 

n=0 

2.2 Cauchy's integral formula 

The definition of operator function by power series can be generalized for wider class of 
functions. To realize this generalization, we use Cauchy's integral formula instead of power 



For example, the exponential 
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series. Cauchy's integral formula states that 

/(„) = -L{!Wi, (2) 

2m J T z — z 

where the integral is a contour integral along the contour T enclosing the point Zo- 

We can define the algebraic isomorphism between an operator algebra and some functions 
[34] . The function f(z) = z corresponds to the operator A. The function f(z — z ) = 
[z — zq)' 1 corresponds to the resolvent operator R(z, A) = (zl — A)^ 1 . If \z\ > r^, where r& 
is a spectral radius: 

r A = lim K/\\A n \\, 

then the resolvent operator exists. The function of linear bounded operator is defined by 

f(A) = ^- I f(z)R(z,A)dz, (3) 
Am J T 

where 

R(z,A) = (zl - A)-\ zep(A). 

Here T = dG G &{A), where cr(A) is a spectrum of operator A, and p{A) C G. For example, 
we can define the operator 

e At = V — A n = — : <b e zt R(z, A)dz (4) 



n\ 2m' 

n=0 



c)G 



that is corresponded to the function exp(zt), and cr(A) C G. 

As the second example, the operator E z that corresponds to the Heaviside function 
9(z — ^o)? where 9(z — z ) = for z > z, and 9{z — Zq) = 1 for z Q < z is defined by 

E z = E{z, A) = ^-id{z- z )R(z , A)dz , (5) 



2ni j r 

and is called the spectral operator. The operator E z can be denoted by 8(zl — A). 

2.3 Spectral representation of self-adjoint unbounded operator 

It is known that spectral function E z exists for all self-adjoint operators A, and 



Ax = J zdE z x, (6) 

where 



\z\ 2 d(E z x 1 x) < oo. 

-oo 



Then, the function f(A) of self-adjoint operator A can be defined by the equation 

r+oo 

f(A)x= / f(z)dE z x. (7) 



As operator A, we can consider the self-adjoint derivative —id/dx. For the father information 
about this approach, we can use Refs. 



3 Fractional Derivatives by Taylor Series 
3.1 Definition of Taylor series 

A one-dimensional Taylor series, which is an expansion of a real- valued function f(x) about 
a point x = a, is given by 

oo 

f(x) = J2fn(x-a) n , (8) 

n=0 

where 

fn = -J (n \a) ) (9) 
n\ 

and f( n \a) is the nth derivative of f(x) evaluated at the point x = a. 

Suppose the function f(x) has all derivatives in the interval \x— a\ < ao, and the condition 

lim i_^Z( a; _ a )™ = o 

ra^oo TV. 



is satisfied, then the series 



/(*) = £— r-(x-a) n (io) 



n=0 



converges to the function f(x) for all intervals [a; — a| < a', where a' < a$. This representation 
of the functions can be used to define fractional power of operator. 

3.2 Taylor series for fractional power of coordinate 

The Taylor series for fractional power of coordinate f(x) = x a about a point x = a > is 

oo 

x a = Y,fn(x-a) n , (11) 

n=0 

where 



n\ ma 



\ n n-a 
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B(a,n) — a(a — l)(ce — 2)... (a — n + 1). 

If m — 1 < a < m, then 

B(a n) = ( iy~m ( a + l - m )m = T ( n - «) r (« + 1) 



(m — a) n _ m r(m — a)r(a + 1 — m) ' 

where 

(z) m = z(z + l)...(,z + m - 1). 

For < a < 1, 

2J(a,n) = {-VT x a{l - e*)n-a = (-l)^^"^. (12) 



Using 

^ S( -i)-(>v- = 5^W> (13) 

"^^C^M,^, (14) 

(_l)*+i a a-*r(n-a) 
(ra-fc)!fc! r(l-a)' 

Equation (I14p represents the fractional power of coordinate as a series of integer powers. 
This representation allows us to define the fractional power of operator as a series of integer 
powers of operator. 

It is known that self-adjoint operators have the real eigenvalues. Using Eq. fjX4j) with 
a > 0, we can define the fractional power of the self-adjoint operator A by 



n 

[x — a t 

we rewrite Eq. (II II) in the form 

oo n 

X 

n=0 fc=0 

where 

C(n, a, a 



n=0 



A« = Y ?^A\ (15) 
^— ' n\a n - a 



For the operator A = —id/dx, we have 

. d \ a B(a, n) ( d 



da;/ ^ n\a n - a \ dx I 

n=0 



or, in the equivalent form 

a oo n 



/ n i — n 



As a result, we get that fractional derivative is defined as a series of integer powers of 
self-adjoint derivative operator. 
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3.3 Examples of computation of fractional derivatives 

Let us consider the fractional derivative (TT6l) of constant c: 



Using 



we get 



If < a < 1, then 



d \ v-^ A(a, n) ( d \ 
t£r / ^ n\a n ~ a V fix / 

7 n=0 v 



-i-^--a] c=(-a) n c, (19) 
ax 



7 n=0 



d\ a „, y> r(n-a) 



' dx) ° ^ I III 



Let us consider the fractional derivative of a power x m . From (fTTj) . 



im = EEH)^ fc > a ' fl )( im ) w - ( 22 ) 

n=0 fc=0 



Using 

= m ( m _ i)...( m _ jfc + i) x m ~ k = m \^ x m - k (23) 

(m — k)\ 

for k < m, and (x m )W = for k > m, we obtain 

7 n=0 fc=0 V ' 

4 Fractional Derivatives by Fourier Series 
4.1 Fourier series 

Fourier series of a function f(x) G Z^f— /,/] is an expansion in terms of an infinite sum of 
sines and cosines. Since sines and cosines form a complete orthogonal system over [—/,/], 
the Fourier series is given by 

oo 

a r {imx\ . . fimx\~\ 

f( x ) = y + [ a ™ cos ) + hn sm )\ ' ( 25 ) 

~ n=l 
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where 

1 f /nnx\ 1 /ixnx\ 

a n = - J f(x) cos — j dx, b n = - J f(x) sin [— — J dx, 

and n is a positive integer number. 

Let us consider the Fourier series for f(x) = \x\ a G L 2 [— /;/] , where a is a positive 
fractional power. The Fourier series of this function for x 6 [— I, +1] is 



where 




The cosine can be represented as a power series. Therefore, equation allows us to present 
the fractional power as a series of integer powers series. Then, we can define fractional 
derivative on the interval [— /, I] as a fractional power of derivative. 



4.2 Fractional derivative 

The fractional power of derivative operator for the interval [— I; I], we can define by 

. d \ a ao / .irk d \ 

=Y + ^ akCOS {-'Td- x l (27) 

7 k=l v 7 

It is known that we can define exp(v4) for the self-adjoint operator A = —id/dx by 



oo 



1 .„ (-i) n ( d 



Using 



-4d/* = Y* A n = Y\-y_tu\ (28) 

^ n\ ^ n\ \dx I 

n=0 n=0 v 



^ = 5(^+0 = 5 E^+E 



n! 14 — ' n! 

n=0 n=0 



we have 



00 00 , , n 00 

y _J_ [(A)" + (-A) n \ = Y ; ,/ A n = Y ^tt^ 2 ™, 

^2(n!) LV ; v ;J ^ 2(n!) ^ (2m)! 

ra=0 v 7 n=0 v ' m=0 K 1 

2)ii 



d \ 1 / c/ 

( US ' _Z / ^ (2m)! V^a; 
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Equation (1671) and (1291) allows us to define the fractional power of operator as a series of 
integer powers of operator. Using (1291) . we rewrite (l27j) as 



(, \ a oo oo / 7 \ 2m / , \ 2m 

% dx) ~ 2 + ^ afc ^- (2m)! I I ) \dx 
' k=l m=0 v 1 v 7 v 



or by the equivalent equation 



where 



Let us compute the coefficients a k by 



1 



ak = I y a cos(Tcky)dy. 

As a result, we obtain 

_!_ Q / 2- a n- 1 / 2+a k a (a + 3) sin(7rA;) 



a k 



2 a n~ 1/2 ~ a r a k- 



(a + !)(« + 3) 



+ 2- Q 7r- 1 / 2+a A; a [7rfccos(7r/c) - sin(Trfc)] + 2-°VkaL(a + 1/2, 3/2, ttA;) sin(7r£;) 



(30) 



7 \ a oo / 7 \ 2m 

/ m =0 v 7 



oo / , 2m 



a + 1 a + 1 

2- a [vrA;cos(7rA;) - sin(7rA;)]L(a + 3/2, l/2,irk) 

where ^, -z) is the Lommel function |38j . 
4.3 Complex Fourier series 

The real- valued function f(x), which is defined on [-L/2,L/2], can be presented by 

oo 

/(*) = E /ne i(27rn/i)s , (33) 

n=— oo 

where 



f n = Y f{x)e-^ nx ' L ^dx. (34) 

^ -/-L/2 
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The operator function f(A) is defined by 

oo 

f(A) = fne i(2 ™ /L)A , (35) 

n=— oo 

where f n are defined in (134"1) . and 

= £ i k ( - 27rn { L)k A k . (36) 

fe=0 

Substitution of ( 1361) into ( 1351) gives 

n=— oo fc=0 

Then the fractional power of self-adjoint derivative on [-L/2, +L/2] can be presented by 



hs) =EE/(«.*)(|)' 

v y n=-oo fc=0 v 



(3* 



= +1 / |x| Q cos(nx)cte. (39) 

Equation (|38|) represents the fractional power as a series of integer powers. This representa- 
tion defines the fractional power of derivative as a series with integer powers of derivatives. 



where 

(2vm^ r- L ^ 



5 Fourier Transform and Fourier Integral 

5.1 Fractional derivative by Fourier integral 

Let us consider a function f(x) with n variables x. Suppose Ai, A 2 ,...,A n are n elements of 
commutative operator algebra A. For example, A is an algebra of operators in linear space. 
We denote by / the Fourier transform for f(x): 



f(y) = I f(x)e lxy dx, xy = x lVl + ... + x n y n , (40) 

where 



m = J^J Jiu)< ,U!U ■-■ !h r '<l <J = ^ / dyf{y)jy A . (41) 
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The operator function f(A) of elements A\, A n is defined by 

f(A) = J^J f(yy^ + - + y^dy = — J dyf(y)^ A . (42) 
Substitution of Eq. flUD into Eq. (020 gives 

-1 /*+oo i /»+00 /»+00 

/(A) = — — / f(y)e^ A dy = — — / dW dxf(x)^ A ' . (43) 

In general, we must have the exact definition of these integrals and description of possible 
class of symbols and algebras of elements A±, .., A n . 

To define the fractional power of the operator A, we use (142]) in the form 

-i /*+oo /*+oo 

A Q = — — / dx dy\x\ a e iy{A - x) . (44) 
(27t)«7_ 0O yM { } 

This equation can be considered as a definition of fractional power of operator A. 
For f(x) = \x\ a , where a ^ —1, —3, we have |39j : 

/+oo 
dx\x\ a e~ ixy = -2sin(7ra/2)r(a + l)|y|" a_1 . 
-oo 

Then Eq. (jUj) gives 

= _2^w2)r (a+ i) r- dyW «-^ (45) 

(2tt)" 9191 ^ ' 

For the self-adjoint derivative operators 

, . d .3 
A\ — — 2— — j ••> A n — —%- 



OX\ ox n 
equation fj4"5l) is 

As a result, we obtain the fractional derivative operator as a fractional power of self-adjoint 
derivative operator. 
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5.2 Fractional power of self-adjoint derivative operator 

Let us consider the self-adjoint derivative operators 



d 

D x = — i— 
ox 



d 







—i- 



dx\ ' ' dx n 



(47) 



It is easy to prove that f(p) = p a 6 S^IR™). Here, 5*°° is the space of symbols that are 
slowly growth on the infinity This space is defined as 



u z n k s?(R n ), 



S°°(R n ) 

where S]° (M n ) is a space of functions from the class C k (R n ) with the norm 

ll/ll S ^) = sup(l + |x| 2 )' /2 (Ell/ (a) ^)ll)- 

1 ran V ' / 



|a| 



The space 5°°(M n ) is defined by 



3r, Vs 



sup{ (i + \x\y 



} < oo. 



The fractional powers of the operator (|4Tj) are elements of algebra /^(S 100 ,^ 00 ) of all 
continuous linear maps of the space S' 00 (R' 1 ). Then, the fractional derivative operator 



Di 



.d_ 

dx 



(4£ 



is a fractional power of self-adjoint derivative operator. The operator (|48j) acts on the 
arbitrary function u(x) G C°°(lR n ) by 



— %- 



d_ 

dx 



u(x) = F p ^ x p a Fy^ p u(y), 



where 

Fy-> P u(y) 

is the direct Fourier transform, and 



1 

2tH 



re/2 



n/2 



e ipy u(y)dy, 



(49) 



(50) 



(51) 



is the Fourier transform. 
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Proposition. The operator 

D °=(-Cf (52) 

has the symbol 

symb{D«}(p) =p a . (53) 
It is not hard to prove this proposition. Using Eq. (I49p . we get 



D a x e p {x) = [ -i-^Je p (x) = ) e*« = F^F^.e*» 

= F Z ^ X (z a {2iri) n / 2 6{z-p)) = 
= F z _ x (p a (2m)^ 2 5(z - p)) = p a e ipx = p a e p (x). 
As a result, we obtain 

D a x e p (x)=p a e p (x). (54) 
Multiplying both sides of (1541) on e_ p (x), we get e_ p (x)D x e p (x) = p a that proves (1531) . 

6 Fractional Derivative by Quantization Map 

6.1 Quantization procedure for coordinate representation 

Let us consider the quantum mechanics in coordinate representation. It is known that 
quantization Q is a linear map of coordinate q and momentum p into self-adjoint operators 

Q(q) = q = q, Q(p) = p = 0(i) = /. (55) 



Using linearity of quantization map [50], SJJ, we get 

Q(aq + bp) = aq + bp. 

Obviously, we have 

Q([aq + bp] n ) = [aq + bp] n . 

Using the power series 

oo ^ 

exp(x) = V -x n , (56) 

z — ' TV. 

n=0 

we get 

Q\exj>[(i/h)(aq + bp)]j — exp[(i/h)(aq + bp)). (57) 
This allows us to define the function of operators q and p by using the Fourier transforms 

HQl SH S21 sai sn SHI 
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6.2 Weyl quantization 

Canonical quantization defines a map of real-valued functions into self-adjoint operators 
[US HU H21 H31 HH 113 06] . A classical observable is described by some real- valued function 
A(q,p) from a function space M.. Quantization of this function leads to self-adjoint operator 
A(q,p) from some operator space M.. 

Let us consider main points of the usual method of canonical quantization [101 |4"2"| H5"| |4"6] . 
Suppose qk are canonical coordinates and pk are canonical momenta, where k = 1, n. The 
basis of the space M. of functions A(q,p) is defined by functions 



W(a, b, q, p) = e (^X<^) , aq = J2 Wk • (58) 



fc=i 



Quantization transforms coordinates and momenta pk to operators and pk- Weyl 
quantization of the basis functions (1581) leads to the Weyl operators 

n 

Q (W(a, b, q,pfj = W(a, b, q, p) = e ^ a ^ , aq = ^ a k q k . (59) 

k=l 

Operators (1591 form a basis of the operator space Xi. Classical observable, which is charac- 
terized by the function A(q,p), can be represented in the form 

i r+oo 

A(q,p) = J A(a,b)W(a,b,q,p)d n ad n b , (60) 

where 

i r+oo 

A( a ,b) = —— A{a,b)W{a,b,q,p)d n ad n b , (61) 



(27TH) n 

i.e., A(a,b) is the Fourier image of the function A(q,p). Quantum observable A(q,p), which 
corresponds to A(q,p), is 

Q (A(q, p)) = A(q, p) = £°° A(a, b)W(a, 6, q, p)d n ad n b . (62) 

This formula can be considered as an operator expansion for A(q,p) in the operator basis 
(155)1 . Substitution of ([EE]) into ((52} gives 

1 /»+oo /*+oo 

^>P) = 7T1^ / / d n qd n pA(q,p)W(a,b,q-qi,p-pi). (63) 

The function A(g,p) is called the Weyl symbol of the operator A(q,p). Canonical quanti- 
zation defined by (I6"3"|) is called the Weyl quantization. The Weyl operator (1591) in formula 
( 163]) leads to the Weyl quantization. Another basis operator leads to different quantization 
scheme [42J. 
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6.3 Fractional derivative by Weyl quantization map 

Let us consider a quantization map of real-valued function f{p) = \p\ a into self-adjoint 
operator. Quantization of this function leads to some self-adjoint operator f(p) from the 
operator space Ai p , where p^ = —id/dxk and k — 1, ...,n. Then 

D: = f(p)=p a = (-zd/dx k r. 

The basis of the space M. p is defined by functions 



W(a,p) = e**, ap = J2 a kPk- (64) 

fe=i 

Quantization maps pk into pk = —id/dxk- Weyl quantization of the functions (|64|) leads to 

n 

Q (W(a, 6)) = W(a, p) = e iap , ap = J2 a*P*- (65) 



fc=i 



The operators fl65|) form a basis of the operator space Ai p . Using Fourier transform, the 
function f{p) = \p\ a can be presented by 



fiP) = J^y^J f(a)W(a,p)d n a, (66) 

where 

1 P + OO 

/(«) = TTTT^ / /(pMM a>, (67) 



(2tt)™/ 2 

i.e., /(a) is the Fourier image of the function f(p). Quantum observable f(p), which corre- 
sponds to f(p), is defined by the formula 

Q(f(pj) = f(P) = j^f 2 f(a)W(a,p) d n a. (68) 

This formula can be considered as an operator expansion for f(p) in the operator basis (159j) . 
From fl67|) and fl68|) . we obtain 

i /*+oo /*+oo 

/(P) = 7T1^ / rf " a / dnp f(p)W{a,p- P I). (69) 

</-oo J —oo 

The function /(p) is the Weyl symbol of the operator f(p). 

For f(p) = \p\ a , where a; is a positive real number, we obtain 

= wr L d " a L d " p mVia - p - piy (70) 

As a result, we have the definition of fractional derivatives on W 1 as a fractional power of 
self-adjoint derivative. 
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7 Fractional Stability 



In this section, we use the fractional generalization of variations of variables. Fractional 
integrals and derivatives are used for stability problems [HI HH [5(S [51] . In this paper, we 
consider the properties of dynamical systems with respect to fractional variations [25]. We 
formulate stability with respect to motion changes at fractional changes of variables. Some 
systems can be unstable "in sense of Lyapunov", and be stable with respect to fractional 
variations. 



7.1 Fractional variation derivative 

Let us consider dynamical system that is defined by the ordinary differential equations. 
Suppose that the motion of dynamical system is described by the equations 

^y k = F k (y), k = l,...,n. (71) 

Here yi, ...,y n be real variables that define the state of dynamical system. 

Let us consider the variation 5y k of variables y k . The unperturbed motion is satisfied to 
zero value of variation 5y k = 0. The variation describes that as function f(y) changes at 
changes of argument y. The first variation describes changes of function with respect to the 
first power of changes of y: 

8f(y) = D 1 y f(y)dy, (72) 

where 

yJ \yJ Qy 

The second variation describes changes of function with respect to the second power of 
changes of y: 

5 2 f{y)=Dlf{y){dyf. (73) 

The variation 5 n of integer order n is defined by the derivative of integer order Dyf(y) = 
d n f/dy n . 

Let us define the variation of fractional order as a fractional exterior derivative of the 
function (zero-form) by the equation 

S*f = D a y f {8y) a , (74) 

where Dy is a fractional derivative with respect to y. 

The fractional variation of order a describes the function f(y) changes with respect to 
fractional power of variable y changes. The variation of fractional order is defined by the 
derivative of fractional order. 
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7.2 Equations for fractional variations 

Let us derive the equations for fractional variations S a y k . We consider the fractional variation 
of equation flTTj) in the form: 

6 a j t y k = 8 a F k (y), k = l,...,n. (75) 

Using the definition of fractional variation ( 1741) . we have 

5 a F k (y) = [D«F k ](5 yi ) a , k = l,...,n. (76) 
From Eq. ( |76l) . and the property of variation 

sa j t ^jf^ < 77 > 

where y k = y k {t, a), we obtain 



d 
It 



-8 a y k = [D«F k ]{8 yi r, k = l,...,n. (7* 



Note that in the left hand side of Eq. ( 1751) . we have fractional variation of S a y k , and in the 
right hand side - fractional power of variation (5y k ) a . 

Let us consider the fractional variation of the variable y k . Using Eq. (J7JJ, we get 

S a y k = [D«y k ] (5 yi ) a , k = l,...,n. (79) 

For the Riemann-Liouville fractional derivative, 

D a m y k * 

if k ^ I. Therefore for simplification of our transformations, we use the fractional derivative 
CIS db fractional power of derivative 

D a yi y k = b kl D a yi y u (80) 

where S k i is a Kronecker symbol. Substituting Eq. (!80l) into Eq. (1751) . we can express the 
fractional power of variations (Sy k ) a through the fractional variation 5 a y k : 

{5y k ) a = Dly k 5 a y k . (81) 

Substitution of Eq. (ED) into Eq. flZSJ) gives 

jfy k = Ktt] [L%F h ] 6"y, (82) 
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Here we mean the sum on the repeated index I from 1 to n. Equation (1821) is equations for 
fractional variations. Let us denote x k the fractional variations S a y k : 

x k = 6 a y k = [D a yk y k ] (5y k ) a . (83) 

As a result, we obtain the differential equation for fractional variations 

^x k = a k i(a)x h (84) 

where 

o«(a) = [D» B ] D«F k . (85) 

Using the matrix X 1 = (xx, x n ), and A a = \\a k i(a)\\, we can rewrite Eq. (1841) in the 
matrix form 

Jt x - KX - (86) 

Equation flHBT) is a linear differential equation. To define the stability with respect to frac- 
tional variations, we consider the characteristic equation 

Det[A a - XE] = (87) 

with respect to A. If the real part Re[X k ] of all eigenvalues A for the matrix A a are negative, 
then the unperturbed motion is asymptotically stable with respect to fractional variations. 
If the real part i?e[Afc] of one of the eigenvalues A of the matrix A a is positive, then the 
unperturbed motion is unstable with respect to fractional variations. 

A system is said to be stable with respect to fractional variations if for every e, there is 
a S such that: 

\\S a v(to)\\ < => \\S a y(a,t)\\ < e Vt G R+. (88) 

The dynamical system is said to be asymptotically stable with respect to fractional variations 
5 a y(t, a) if as 

t^oo, \\S a y(t,a)\\ ^0. (89) 

The concept of stability with respect to fractional variations is wider than the usual 
Lyapunov or asymptotic stability. Fractional stability includes concept of " integer" stability 
as a special case a — 1. Some systems can be unstable with respect to first variation of 
states, and be stable with respect to fractional variation. Therefore fractional derivatives 
expand our possibility to research the properties of dynamical systems. 
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